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Solution of string equations for asymmetric potentials 


Patrick Waters 1 ’* 


Abstract 

We consider the large N expansion of the partition function for the Hermitian 
one-matrix model. It is well known that the coefficients of this expansion 
are generating functions F ^ for a certain kind of graph embedded in a 
Riemann surface. Other authors have made a simplifying assumption that 
the potential V is an even function. We present a method for computing F^ 
in the case that V is not an even function. Our method is based on the string 
equations, and yields “valence independent” formulas which do not depend 
explicitly on the potential. We introduce a family of differential operators, 
the “string polynomials”, which make clear the valence independent nature 
of the string equations. 
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1. Introduction 

We consider a partition function defined by 


Z n ,N{t) = / exp 


-nJ2v( K;t) 


i= 1 


n (A, — A,) 2 d"A. 


We assume that the potential function V is a polynomial: 

d 


r(A;f)=-A 2 + 

3 =1 
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For the integral in ([!]) to converge, it is necessary that the polynomial de¬ 
gree d is even, and that the parameter td is such that the real part of td is 
positive. By replacing the domain of integration R n with T n , where T is an 
appropriately chosen contour in the complex plane, it is possible to remove 
these restrictions [3]. 

The partition function in ([!]) is the normalization constant for the induced 
measure on eigenvalues arising from the following probability measure on 
n x n Hermitian matrices: 


dP n (M ) = —?:-exp [~NV(M- 1)] d M. (3) 

Z n ,N\t) 

The factor d M is Lebesgue product measure n 2 real degrees of freedom of 
M. 

We are interested in the partition function (JT[) because it has a topological 
expansion whose coefficients are generating functions for a certain kind of 
embedded graph call a map. Letting n, N —» oo with the ratio n/N = x held 
fixed, we have 


log 


Z n ,iv(t) 
Zti,n( 0) 


J2 F{h) (x,t)N 2 ~ 29 , 

9 > 0 


p(s) 


X 


Faces (r) 


re 


genus g maps 
vert’s of valence < d 


n 

3 = 1 


-tj) 


FT) 


B(r)i 


( 4 ) 
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Here V^ (r) is the number of j-valent vertices of T: that is, vertices at which 
exactly j edges meet. We use ~ to mean that for any g > 0, truncating 
the series after terms of order N 2 ~ 29 results in an approximation with error 
0(N l ~ 29 ) as n, N —y oo with n/N = x. 

A map of genus g is an equivalence class of labeled graphs embedded 
in genus g compact connected oriented surface. Two embedded graphs are 
equivalent if an orientation preserving homeomorphism of the surface takes 
the vertices and edges of one to those of the other, preserving labels. A map 
is equipped with the following labels: each of the map’s kj vertices of each 
valence j has a unique label 1 ,,kj] and also for each vertex, one incident 
edge is marked with an arrow pointing away from that vertex. This marking 
of edges can be represented as a function / from the vertex set to the edge 
set such that each vertex V is incident f(V). 
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The expansion Q was used by Bessis, Itzykson and Zuber [3]. Ercolani 
and McLaughlin [9] gave a rigorous proof that an expansion of the form 
Q exists, and proved analytical properties of the expansion. In particular, 
expansions of derivatives of log Z n)N (t) / Z n ^ N {{)) with respect to g can be 
computed by termwise differentiation of the series Q. The same is true 
for differentiation with respect to x, but for more complicated reasons; see 
section 4.1.3 of na for a proof that this property is equivalent to the string 
equations. 

This paper addresses the problem of computing the generating functions 
F^, a problem which has been studied by several authors jH} EJ IU 113; HI 
ei eu sun mg. These articles give a variety of methods for computing F ^ 9 \ 
and yield formulas given in terms of a variety of choices of variables. In this 
article we focus on valence independent formulas. Formulas of this kind allow 
derivatives with respect to x but forbid explicit dependence on the potential 
V or its parameters tj. For example we find that 


F m = - 4 log i log ((' u') 2 z - (o') 2 ) , 


( 6 ) 


where z and u are defined in equations (16 17) and the notation ' indicates 
d x . Other authors assume that the potential V is an even function; this 
assumption significantly reduces the difficulty of computing F^ 9 \ for example 
with this assumption we have u = 0. Thus our formula (|6]) is new and extends 
formulas given by other authors to the case of an asymmetric potential. The 
valence independent approach was first used by Shirokura hi. 

For the purpose of computing the generating functions F^ 9 \ it is difficult 
to work with the partition function directly. Szego’s relation (J7]) and the 
Hirota-Szego relation ([8]) give two fundamentally different relations between 
the partition function and recurrence coefficients for a family of orthogonal 
polynomials: 


J n -\-1 


,N Z n —l 


Z'Ln 


r n,N 


n + 1 

^n,N • 

n 


i d 2 


log 


Z n ,N(t ) 
Z h ,n{ 0) 


(7) 

( 8 ) 


Strategies for computing F^ may involve using string equations to compute 
higher order asymptotics of the recurrence coefficients r n> jv, then lifting this 
information to the level of the partition function Z n ^. For example [3] 
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uses 0. pm uses (J8]) and [T] uses Blcher-Its deformation. We will use 0; 
the mathematical rigor of our use of this formula was established in section 
4.1.4 of [15]. Strategies for computing F ^ which are based on loop equations 
M avoid the recurrence coefficients r n , but involve computing higher order 
asymptotics of random matrix correlators. It is also possible to compute F 
directly from the Riemann-Hilbert problem for orthogonal polynomials [T5] . 

Let us now define the coefficients r n and explain the string equations 
which govern them. The monic orthogonal polynomials p n ^N arise from the 
following (non-Hermitean) inner product and orthogonality relation: 


(F(A),G(A)) = [ F(X)G(X) exp (-NV (A; t)) dX. (9) 

Jr 

(Pu,n{ A), X m ) =0 for all rn < n. (10) 


Thus the polynomials p n can be constructed by the Gram-Schmidt process. 
As is well known [Tjj, the polynomials satisfy a three term recurrence: 


Apn,w(A) — Pn+l,v(A) + S ni JvPn,v(A) + ^n,A r Pn-l,Af(A). (11) 


This defines the recurrence coefficients r n . The three term recurrence 
can be expressed in terms of a Jacobi matrix L: 


(ID 


/ s 0 1 0 0 \ 

ri si 1 0 

0 r 2 s 2 1 


' 7 


( 12 ) 


Let p = (p n (A))n>o be the column vector of all orthogonal polynomials for 
the potential V ; then the three term recurrence can be encoded as a matrix 
equation 


A p =Lp. (13) 

Following Blccher and Deafio [4], we use the string equations 

0 =V\L) n ^ (14) 

£ =V\L) n , n _ i. (15) 
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A differenced form of these equations was used in ra- Equation ( |14[ ) can be 
derived as follows. Integrating by parts we find that 


0 = J p' n (X)p n (X) e~ NVW dX 

=N J p n (X) 2 V\X) e ~ NVW dX 
=N\\p n \\ 2 V'(L) ntn . 


An analogous calculation, but starting with n \\p n -i IT = f p' n (X)p n - i(A)e NV W dX 
establishes (15). 


We now explain how the continuum limit of string equations can be used 
to compute the generating functions F^ 9 \ The following asymptotic expan¬ 
sions for recurrence coefficients with shifted indices was proven in [10] : 


r n +k,N ~ exp (kN 1 d x ) ^ z g (x, t)N 2ff , 


g> o 


s n +k,N ~ exp (kN 1 d x ) ^ u g (x, t)N 

g> o 


-g 


(16) 

(17) 


Thus we have defined z g and u g to be the asymptotic coefficients in the 
expansions (16 17). We will use the abbreviations z = z 0 and u = Mo- 


Inserting these expansion in the string equations gives the continuum limit of 
the string equations. By extracting terms at order N~ 2g from the continuum 
limit of string equations, one can solve for z g and u^ g in lower order terms. 
Inserting the expansions (16 17) and the topological expansion Q in the 
Szego relation ([T), one may solve for F^ in terms of the variables z g ', g' < g: 


F {g) _ (1 - 2m)B 2m Q 2 m -2 


m =0 


(2m)! 


Zg—m 


(18) 


The constants i? 2 m are Bernoulli numbers, and the cumulants z g are defined 
by 


z g N~ 29 ~ log x-'zgN- 29 . 
g> o 9>o 
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In 0, formula ( |18| ) was obtained using the Eulcr-MacLaurin formula. As was 
noted in [I], further justification is required to make this approach rigorous. 


A rigorous proof of (18) has been given in section 4.2.1 of [15]. 

The rest of our paper is organized as follows. In section 2 we prove a 
theorem concerning the structure of the string equations; theorem [T| describes 
the terms arising in the string equations as a sum of differential operators 
(which we call string polynomials) acting on a residue. The reason that our 
method results in valence independent formulas is essentially that the string 
polynomials do not depend on the potential. In section 3 we demonstrate 
our method by computing F^. In section 4 we prove a structure theorem for 
valence independent formulas for u g and z g . To compute using formula 
(18) it is necessary to evaluate an indefinite integral with respect to x. It 


turns out that in the cases <7 = 1,2 this integral can be computed in closed 
form, but we have no proof that this happens for all g > 0. 


2. String polynomials 


Let us sketch an outline of this section. In 2.1 we state theorem [lj assert¬ 


ing that the ^-derivative structure of the string equations can be expressed 
in terms of a family of valence independent “string polynomials”. We then 
reduce theorem [T| to lemma [lj which essentially allows us to work with each 
time variable separately. In sections 242 - 2 A we prove lemma [l] Our proof 


gives a formula for the string polynomials, which is stated in lemma [3j 
2.1. String polynomials 

In this section we will use the following notations 

r n ,N = y ^z g N~ 29 , 

h> 0 

■V.'Y ~s = u 3 N~ g . 
g> 0 


( 21 ) 

( 22 ) 


That is, we are defining r and s to be the formal series for the asymptotic 
expansions of r n ^ and s n Note that ~ has the same meaning as in (j5|; 
let us give two examples marking the distinction between ~ and =. First 
exp(-AQ ~ £,>o 0 x N 9 = 0, where the sum converges to the “wrong” 
limit; and second J 0 °°(l + r/A^) _1 exp(—r) dr ~ ^> 0 (—1 ) 9 g\N~ 9 , where the 
sum diverges (this example is the Sticltjes function, see [3]). 
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We will use the notation y/r for the asymptotic expansion of y/rf, where 
we choose the branch of the square root that is positive when the potential 
is real. The square root is natural because if one works with orthonormal 
polynomials P n (with positive leading coefficients) instead of monic orthogo¬ 
nal polynomials p n , the three term recurrence becomes A P n = y/r n+ \P n +i + 
s n P n + y/PnPn-i- Here r n and s n are the same coefficients as in the case of 
monic polynomials. 

The variables r and s allow us to analyze the ^-derivative structure of 
the string equations before stratifying terms by genus. We will use integer 
partitions as multi-indices to describe differentiation: 


len(A) 

a x j(x) • II dtf(x). (23) 

i =1 


Furthermore, since x-derivatives in (16 17) are always accompanied by pow¬ 
ers of TV -1 it natural to differentiate with respect to n — xN, that is 


d n f(x) := -/V _1 /'(x). Thus for example <9„ +2 S' = N~ 3 S'S”. 

We will use the notation [h p ] for the h p coefficient of a Laurent series in 
powers of h. That is, 


m)^{ G w h -’ f. 

In section \2A\ we will prove the following theorem. 

Theorem 1 (Existence of string polynomials). There are polynomials (with 
coefficients in the ring Q[r]j P^ (d r , d s ) indexed by integer partitions X,rj 
such that the continuum string hierarchy can be expressed as 

o ~ d " r (&> d r ) • [h°]V'(h + s + rh~ l ) (24) 

X ~ ^2 &n s dy d r ) ■ [h°]V'(h + s + rh~ l ) (25) 

A,?7 

The string polynomials Py rj ° r ^ do not depend on the potential V. Further¬ 
more, they can be chosen to have polynomial degree at most one in r, and 
are uniquely determined by this condition. 
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We now reduce theorem Q]to a lemma which will allow us to work with one 


power of the Jacobi operator L at a time. Define the polynomials P^ j5 
by 


P 


(&) 


A ,T],J 


(■ L ' J X )n,n ~ Y Piij (s, r) d*s d*r, 


\,T) 


{L j 1 )„,„_! ~ Y P \l,j ( S > r ) d n s d n r - 


(26) 

(27) 


Unlike the polynomials P^ or b \ it is clear that the polynomials 
For example, if we let J — 3, then (26) becomes 


(a or b) 


exist. 


p n,n ~ r n +1 + 4 + r « 

~ (s 2 + 2r) + iV _1 (2ss' + r') 

+ N- 2 ({s') 2 + ss"+ l -A+.... (28) 


The hrst several polynomials P\^A s i r ) can be immediately read off from 
the above example. For example P^ 3 (.s, r) = 2s, and P[°+] (j) 3 (s, r) = 1. 

Lemma 1. For all \,r] there exist polynomials P ^, P P q , such that the fol¬ 
lowing identity holds for all J: 

P \lA s i r ) = p x$( d s, d r) ■ [h°](h + s + rh~y~\ (29) 

P \lj( s ,r) =Pj®(d 8 ,d r ) ■ [h~ 1 }(h + s + rh- 1 )- 1 - 1 . (30) 

We will prove lemma [T| in section |2.4| 

Proof that lemma [I] implies theorem [Ij First consider the string equation 0 = 





V'(L) njn . We must show that the right hand side is as given in theorem [TJ 
V'(L) n ,n = Y + [1 ^ J= 2])(^ J " 1 )n,n 

J> 1 

^Y 9 n sd n r Y J ^ J + t 1 if J = 2 ])^S,j( S ’ r ) 

A ,?7 J>1 

=E^s>E^+i 1 ifJ = 2 ]) 

A ,77 J>1 

x P^(d s , d r ) ■ [h°](h + s + rh~ l ) J ~ l 

= Y d n sd n rP X^( d ^ d r) ’ [h°]V(h + S + rh~ l ). 

A,)? 

An analogous calculation applies for the other string equation x = W(L) n>n _ 1 . 

d 


2.2. Motzkin path formulation of string equations 

A Motzkin path of length £ is a function p from {0,1, 2,..., £} to Z such 
that p(m + 1) — p{m) e {—1,0; 1} f° r each m. We also think of p as the 
piecewise linear function obtained by joining points on the graph of p with 
line segments. Because the matrix L is tridiagonal, the entries of a power of 
L can be expressed in terms of a sum over Motzkin paths. Let Mq J 0 ^ be 
the set of Motzkin paths of length J — 1 that start and end at height 0 (that 
is, p(0) = p(J — 1) =0). Then 


Ltg= v co) 


(31) 

P ^ M oV 



len(p) —1 I 1 

if p(i + 1) > p(i), 


^i.P) • \ Sn-\-p(i),N 

if p(i + 1) = p(i) . 

(32) 

i=0 1 „ 

f ' n+p(i),N 

if p(i + 1) < p{i) 


We call the quantity C(p) the contribution of the Motzkin path p. The string 

equations can be expressed in Motzkin path form: 


0 + Y^ 

E c(ri, 

(33) 

i 1 

56 M 3 n -f 


n \ 

X^n + Yltj 

Y c (p)- 

(34) 


j P eM n,n-i 
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which defines the polynomials P ^, it is clear that there exist polynomia 
C\. rj {p) such that 


We wish to use the Motzkin path form of the string equations to state a 
formula for the modified string polynomials P[° b \ From equation (26), 


C(p ) =J2 c *M d nsdy 


(35) 


The quantities Cx, v (p ), which are polynomials in r and s, can be easily found 
by substituting the asymptotic expansions for index-shifted recurrence coef¬ 
ficients and collecting terms. Thus 


=£& 


A ,rj 


E 


P £M 0fi 1] 


■ J2 Cx M d n s9 nP 


and therefore we have 


<b= E 


PG< 0 1 


d*b= E C >M- 


(36) 

(37) 


p^K-i 


For example, consider the Motzkin path p shown in figure 1. The path p 
contributes C(p) = s n+ ir n+2 r r)+ ir n s^_ 1 to the quantity Lq 0 , calculated as a 
sum over Motzkin paths according to formula (31). To find the contribution 
C(p] 10,1, 4>) of this path to the string polynomial P 2+1 <*, we substitute the 


series 


s n +k,n + — N V + —N 2 s" + ..., (38) 

r n +k,n ~r + Y\ N ~ lr ' + \\ N ~ 2 ' r " + ''' ( 39 ) 
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Figure 1: For this Motzkin path C(p) = a n+ i&^ +2 6^ +1 6: 


n u n-1 • 


in the formula for C(p ) and read off the N 1 s' term: 

C\p) =s n+1 r n+2 r n+1 r n s 2 n _ 1 

=s 3 r 3 + IV -1 (— s 2 r 3 s' + 3 s 3 r 2 r') + 0(N~ 2 ) 

=k Ci,</),io(p) = - sV 3 . 

Notice that in this example, computing the contribution of a path to a string 
polynomial amounts to selecting one flat step of the Motzkin path from which 
to choose the N~ 1 d x term in expansion (38); from all other steps we choose 
the leading order term s or r. 

We now give a formula for modified string polynomials in terms of Motzkin 
paths. Define a marking of a Motzkin path p with partitions A, p to be a pair 
of functions </>_►, 0^, where: 0_j. assigns to each part of A the left endpoint 
of a distinct flat step of p , and assigns to each part of 7] the left endpoint 
of a distinct downward step of p. The contribution of the Motzkin path p to 
can be computed as a sum over all markings of p with A and p. Also, for 
a path p let —> (p),\ ( p ) respectively be the number of flat and downward 
steps of p. 


Lemma 2. 

>(“) 


pi x a lA s p)= s ~ 

P eM ofo _1) 


•(p) r \(p) 


n 


P(<M0) Ai TTP(K( i ))’ w 


AJs 


n 


Vi- r 


(40) 


Proof. Notice that the dependence of the right hand side on A and r] is 
through the conditions on the markings explained above. We com¬ 

pute the modified string polynomial according to formula (36). For each 
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downward step of the Motzkin path, we must sum over possible choices of 
a term in the expansion r n+p p) = 'Yh m p{i) m d™ r /m\. Similarly for each flat 
step of the path. The sums over the markings (J)\ and </>_> encode these 
markings. □ 


2.3. Example P^ 

In this section we give a formula for the string polynomial P ^ in the 
case that A is the empty partition, and 77 has a single part of size m. We will 
then restrict to the case m — 2 and compute P^l explicitly. Our proof of 
theorem [T| in the next section simply generalizes the calculation here. 

Appealing to formula (40), let us identify the conditions on the sum over 


< t >_> and (f)\. Since the partition A = (j) is empty, there is no sum over choices 
of The function (j)\ assigns the single part m of the partition rj to one 
of the downward steps of the Motzkin path p. This is the unique downward 
step i from which we choose the N~ m d™ term from the expansion on r n+p u\. 
Let ii,j 1 , k\ be the numbers of upward, flat, and downward steps of p before 
this special step, and 12 , J 2 , be the numbers of upward, flat and downward 
steps after it. Clearly we have the following relations 


E 


a =1 


j a + k a J 2 


E 


k n =1 


a=l 


(41) 

(42) 


Notice that from any given starting height, the number of Motzkin paths 
with numbers i,j, k of upward, flat, and downward steps will be given by the 


trinomial coefficient 
calculated as follows: 


* 

i, j, k 


_ (i+j+k)\ 
i\j\k\ 


. Then the string polynomial can be 


E 


s 3i+h r 


kl+k2 (h - h y 


^OL~\~ja.~\~k OL — J— 2 
'y y q, ia k(x = 1 


X 


ml 

* \ / * 
h, ji, h) V*2, 32, k 2 J' 


(43) 


Notice that the height before the downward step from which we take the 
N~ m d™ term is given by the difference between the numbers i\, k\ of upward 
and downward steps. 

We now use combinatorial tricks to evaluate the right hand side in for¬ 
mula (J43|) . Instead of explicitly putting conditions on the sum, we introduce 
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dummy variables £ for the length of the path, and h for the change in height 
of the path. 


il,...,k 2 >0 


(i 1 ~ h) r 


ml 




a=l 


? jot ? k Q 


= [h 1 e J ~ 2 ]\ h ^ h V (?1 ^ m r -v2 


21,...,/C2>0 


= [hH J ~ 2 ]\ h ^ h 


{hid; 


a=1 
2 


hu 


ml 


a= 1 


= [h 1 £ J ~ 2 ]\ h ^ h 


(. hid hl ) r 


mi 


r 1/2 Y\(£\/rh a y°‘(£s) ja (£^/rh a i y 

ii,...,k 2 >0 

1/2 n 

a= 


-r 


1 - £(VrK + S + A/r/r- 1 )" 


^or) jot-) k* 
* 

i'Oti jon K 

(44) 


The order of operations reads from right to left: apply derivatives, set h Q ’s 
equal to h, then extract the series coefficients. Let us comment that in for¬ 
mula (44), the coefficient of £ J ~ 2 must be extracted first, then the coefficient 
of h ] , and not vice-versa. 

For any particular value of m, we are able to evaluate the derivatives 
in equation (44) symbolically, then use familiar Taylor series to extract the 
appropriate coefficients. Consider the example m — 2. Evaluating the two 

(. h 2 - 2 + h~ 2 )£ 2 ^/r 


derivatives in (|44|), we are able to calculate 
>(“) 


pW( s r -J) =\h l £ J ~ 2 1_ (h + h l )£/2 

^’ ’ > [ (1 — £s — £y/r(h + h- 1 )) 3 ' (1 ~ £s — £y/r(h + h -1 )) 4 


+ 


-\hH J - 2 \ ( \d,djf + \sftd,3 ,y - \fi&, 


o -1 


1 - £s - £^r(h + h- 1 ) 
yl'f/r + i ,,l ]!v / r/' + A + x/rh- 1 ) ’-' 

(45) 


= ( jL» 2 + ) [h°](h + s + rh-') J -' 


We have expressed the coefficient Pfj^ j(s,r) as a differential operator acting 
on (h + s + rh~ 1 ) J _1 , exactly as in equation (29). In this case we see that 
the differential operator does not depend on J; theorem [l] asserts that this 
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will be true in general. We have shown that 


P { °\ =-d 2 s +—d r . 

<t>,2 Q l2 r 


( 46 ) 


For the last step in calculation (45), we used the following identities: 


d s [h}-]{y/rh + s + \frh 1 ) J 1 =-d^[h°](Vrh + s + y/rh 1 ) J 1 (47) 

d 2 s [h°}(Vrh + s + y/rh~ 1 ) J ~ 1 (r~ 1/2 d^ + d^j [h°](y/rh + s + y/rh~ 1 ) J ~ 1 

(48) 


These identities are somewhat trivial, but we give proofs in section Appendix 
[A] see lemmas 10 and 11 


2-4- General case 

We now extend the method in the above example to compute P[ a ^(d s , d r ) 
for arbitrary partitions A, rj. Clearly the polynomials can be calculated 

analogously, but replacing Mqq 1 ' ) with Mq J _^ in formula (40). The fol¬ 
lowing lemma gives an explicit (though cumbersome) formula for the string 
polynomials: 

Lemma 3. Let /(A) and l(rj ) be the numbers of parts of the partitions A, 77 . 
Let I be the left ideal of the ring C[r 1 / 2 , r~ 1//2 , d r , <9 S ] generated by rd 2 +d r —d 2 . 
Then P^(d s ,d r ) is the unique element of degree at most 1 in d r within the 
following coset: 


I Pi, (9„ 9 r ) = £ C„( A. r,)Wv(d„ 9^)??* >9« mod I. 


(49) 


The coefficients Wp are defined by 

[£(hd h ) Vq (l - £ {yfrh + s + Vrh -1 )) -1 ] 

<7=0 

fW+Kv) 

= I Vs{d„ d^p) ^ a , (A)+i( , )+1 s/r(h - hrpjlpl) 


(50) 


g^)i-‘('/)(i _ £ 4 - s 4 - yT/r x )) 1 if ^2 v i even 


1 — £ (y/rh + s + a /rh~ l ) 


if ^2 Vi is odd 
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The coefficients C$ are defined by 


IW+Kv) ("V 9-1 hr) — rr V 
J-j- lz^q= 0 ' l <yCJh a CTgJ 

m,o q= 1 


m. 


q- 


= V C„-(A, •]) 


l(\)+l(ri) 

n (w 

< 7=0 


(51) 


The conditions on the sum overm , cr are that the values ofm : { 1 , 2 ,..., /(A) + 
K 7 l)} —t ^ are 7n 1-1 correspondence with the parts of A and rj; and that 
a : { 1 , 2 ,..., Z(A) + l(rj)} —» { 0 , 1 } is sac/i that if a(i) = 0 , then m(i) is a 
part of X; and if a{i) = 1, then m{i) is a part of rj. 


As in section 2.3 we start with the Motzkin path formula (40) for the 


modified string polynomial Pj. To generalize the calculation there, we let 
i a ,ja,k a be the numbers of upward, flat and downward steps between the 


a th and a + 1 st marked steps. To give a formula analogous to equation (43), 


we need to define two auxiliary functions that describe which parts of the 
partitions A, p correspond to which marked steps. Let cr p be the function 
function that takes the value 1 if the p th marked step is a downward step 
and zero otherwise. The function m represents an assignment of parts of the 
partitions A, r) to the marked steps. 


p iiA s ’ r )=Yl S 


IW+Kv) (spq-l ■ , _ y 

yZyo=0 U 0‘J 


m ’ a T, i q~ k q = Kv) q= l 

J2 iq+jq+k q =J-l-l(X)-l{r]) 

iW+Kv) 

- n 

q=0 


m q \ 


; ) (52) 

Jqi roq , 


The conditions on the sum over m, cr are that the values of m : { 1 , 2 ,..., Z(A) + 
l(r))} — > Z are in 1-1 correspondence with the parts of A and r/; and that 
a : {1,2 ,..., Z(A) + l(rj)} — > {0, 1 } is such that if a{%) = 0, then m(i) is a 
part of A; and if cr(i) = 1, then m(i) is a part of r /. Notice that the principal 


way in which the right hand side of formula (52) depends of A, r] is through 
the conditions on the sum over m, cr. 
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By a calculation analogous to the example in section |2.3[ we find that 

1 


P$(s,r) = JJ[h l 


l(A)+l(tj) (sr^q -1 / q _ n V 

n iLtt =0 h a O ha (Tg) 

h a yh <?—1 


m q \ 


x 


n 

<?=0 


1 - l{yjrhq + s + A /rhg 1 ) 


(53) 


The dependence of terms in the above sum on m and a appears only through 
the differential operators acting on the product of generating functions. We 
can thus expand this differential operator as a polynomial in dh x ,..., dh l(x)+l(ri) , 
defining the symbol C#(\,r)) for the coefficients of this polynomial; we have 
chosen to explicitly indicate that the coefficients depend on A ,r]. 


JJ- V,/—7q= 0 n aOh a - &a) 

m,a q= 1 


m. 


<?• 


iW+Kv) 

If tw 

v q =0 


The sum is over tuples v = (vo,V\, ..., vp\) + p v )), but clearly only finitely 
many terms in the sum are nonzero. In terms of the coefficients just defined, 
we find the following formula. 


p( a ) 


(s, r) 


1 

fliv )/2 


E C «-( A > i) 


iW+Kv) 

n 

q =o 


(hd h 


1 — i{^frh + s + y/rh~ l ) 


Following the example in section 2.3, we would like to represent this as a 


differential operator given in terms of d r and d s acting on 1/[1 — f'( v / rh + s + 
-y/rh -1 )]. The purpose of the following lemma is to exchange dh derivatives 


for d s and d s derivatives. 


Lemma 4. For each tuple of nonnegative integers v, then there exists a 
polynomial W${d s , d^) with coefficients in the polynomial ring C[a/t] such 
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that 


a+/3 

n 

q =o 


G hd h y 


1 

1 — £ (y/rh + s + y/rh- 1 ) 


= r m- i (v) Wv ^d s ,d v - r ) 


{ qKX)+ l(jl) _ 1 _ 

1 — £ (y/rh + s + y/rh~ l ) 

B i{X)+m +1 Vr(h~ h^/Kv) 

1 — £ (y/rh + s + y/rh- 1 ) 


(54) 

if Y v i is even 
if Y v i 'is odd 


Proof. First we claim that for any nonnegative integer v there exist rational 
functions R VtP (h), p = 0... v with integer coefficients such that 


(hd h ) 


V _ 

1 



V 


£ 


(U/rTR^K) 

(1 — £ (y/rh + s + y/rh~ 1 ))P +1 


(55) 


In fact furthermore we claim that if v is even, then R V:P is a polynomial 
function of h+h' 1 ; and if v is odd, then R V)P equals h—h _1 times a polynomial 
function of h+h _1 . It is easy to verify this claim by induction on v. It follows 
that 


iW+Kv) 

n w' 


q=0 


1 

1 — £ (y/rh + s + y/rh- 1 ) 



(£B)PR?, p (h) 

(1 — £ (y/rh + s + /rh- 1 ))P +l R)+h , n) ’ 


(56) 


where R$ jP is a polynomial function of h + h 1 is Y v i is even, and h — h 1 
times such a function otherwise. The lemma now follows from because 


f) a f) b _—_ 

s ^1 — £ (y/rh + s + y/rh- 1 ) (1 


£ a+b (h + h- 1 / 

£ (y/rh + s + y/rh~ 1 )) 1+a+b 


(57) 


We note that the asserted divisibility by ds X, ~ l(ri> 1 is because d s is needed 
to increase the exponent in the denominator without increasing the degree 
of the numerator. □ 


It is convenient to express the right hand side of (54) in this form because 
the following identity will allow us to recombine the even and odd Y v i cases 
in the next calculation. 


yfrd s [h p ](h — h /(y/rh + s + y/rh l ) v =p[h p ](yfrh + s + y/rh ^ (58) 
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Formula (58) is a special case of the “integration by parts” lemma proven in 


section Appendix A It follows that 


P^IjM 


= [h\ 




fl(v )/2 


£ c,r‘^w, {a „ m + , + 


v: ^2 Vi even 


iT: v* odd 

= ^Ec iT iv ( ,(s,,s Vf )a? A,+,( ’> ) [ft l| ’ l ](V?/ i + s +v?*- 1 )-'- 1 

V 

= J2 CvWvids, d^)d l s ^d l r ^[h 0 ]{sfrh + s + vAA-y- 1 (59) 

V 

= P^(d s ,d r ) [h°](y/rh + s + a/ r/T 1 ) 7-1 . (60) 


To get (J59[) we applied the “repeated path lowering identity” of section Ap¬ 


pendix A. Clearly, there exists P such that the final expression (60) of our 


formula for string polynomials 
identity” of section Appendix A 


rolds; however by the “derivative swapping 
we can and do choose P to have polynomial 
degree at most one in d r . This completes the proof of lemma [lj and thus 
also theorem [l] 


3. Solution of string equations 


0 : 


The following valence independent formula was known to Bessis et. al. 


x rxi 


= / / log — da:? daq. 

J 0 Jo X 2 


(61) 


We are not aware of any closed form evaluation of the integrals in (61) which 


holds for general asymmetric polynomial potentials. However in the case of 
a potential V(A) = \\ 2 + t\\ where j is allowed to be odd, we have found an 


explicit formula [15]. We consider (61) to be the valence independent formula 


for F'°); thus the real work begins at genus 1. In this section we compute 
using the string equations. 
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3.1. The unwinding identity 

Define two sequences of functions by 




■ = {h Q ]V( m+l \h + u + zh- 1 ) 


h 


flm+l 

u u 

zd™ 


hi \ 

q J [h°]V(h + u + zh~ x ). 


(62) 

(63) 


Notice that the above expressions depend explicitly on the potential V. The 
functions 0 m and 0 m are useful because it turns out that they can actually 
be expressed in a form independent of the potential (“valence independent”) 
are furthmore satisfy the “unwinding identity”. 

Lemma 5 (Valence independence of 0 m ,0 m and the unwinding identity). 
The functions <f > rn , 0 m are determined by the following recurrence 


f 4 > m +1 
\ " 0 ) 71+1 


1 / -zu' 

(z') 2 - z{u') 2 V zz' 



(64) 


along with initial conditions 0o = 0 and ip Q = x. For integers m > 1, we 
have 


Z fm T W 0m 0 m _ 1 

Zu'fpm + z'lfm =Z(j)' m _ 1 


(65) 

( 66 ) 


In the case m. — 1 the right hand sides of the unwinding identities are 
0o = 1 and z4>q = 0. We will use this in section [2] to isolate highest genus 
terms when solving the string equations. The case m — 1 was used in HU, 
and the proof we give below essentially parallels the calculation there; lemma 
[5] extends Ercolani and Pierce’s idea to m > 1. 


Proof. We proceed by induction on m. The base case rri — 0 is the leading 
order equation of the continuum string hierarchy. Let Y — h + u + zh -1 . 


(fr 0n 

(d'x 0r; 


= [/r°]V (m+2) (V) 


1 h - 1 
h 1 


u 

z' 


(67) 


By elementary linear algebraic manipulations, we see that (67) is equivalent 
to the following system; we then apply the path reflecting identity, which is 
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lemma [9] of appendix |Appendix A 


3 x d 
zd x <t>. 


m 

m 


=[h°]V (m+2) (Y) 
= [h°]V {m+2 \Y) 


1 z 1 h 
zhr 1 1 

1 h - 1 
h 1 


zu 


zu 


( 68 ) 


We now combine (67) and (68) to form a matrix equation; once again using 


the path reflecting identity we find that 


^xPm ^x^Pm 
dx'Pm ^^xPn 


U 

z' vl 


[, h°]V {m+2 \Y ) 


1 h 
h z 


(69) 


The hrst column of the above equation is equivalent to the recurrence (64) 


for (p m , ipm- Rearranging equation (69) and expressing the series coefficients 
in terms of <p m and pJm gives the unwinding identity. □ 


3.2. Extraction of generating functions 

So far it has been useful to work in terms of the asymptotic expansions 
s = YPN~ 9u g and r = YPN~ 2a z g . But to solve for z g ,U 2 g we will need to 
extract terms from these expansions. For this purpose, we have the following 
formula. 

Lemma 6. IfG(s,r ) is a polynomial, then 

[N~ p ]G(s, r) — V n^[iV°]3fW’’ ) G(s,>-) (70) 


Proof. Consider the case where G is a monomial s p r q ; we have a contribution 
corresponding to each pair tuples (g a )a=i an d (f/ a )a=i which uniquely deter¬ 
mine a way of choosing a terms u ga N~ 9a from each factor s and z g ' a N~ 29a 
from each factor r. The partitions A, g are the multisets corresponding to 
the tuples g,g'. The number of tuples corresponding to a partition in this 
way can be expressed in terms of factorials in such a way that we arrive at 
the lemma statement. □ 

The sum is over partitions A and g, and by l( A) we mean the number of parts 
of the partition A. Thus we are able to calculate that 

[N~ 1 h°}V'(y/rh + s + y/rh^ 1 ) —Uid s [h°]V' (y/rh + s + y/rh~ l ) 

=UiPi, (71) 
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and also that 


[N- 2 h°]V\V^h + s + V^h- 1 ) 


= yzid r + u 2 d s + 

=-101A + U 2 (j) 1 + ^«i02- 


[h°]V\y/rh + s + y/rh x ) 
(72) 


3.3. Example: solution at order N 2 

Using onr table of string polynomials, we can now write down the equa¬ 
tions at order N~ 2 from the continuum string hierarchy. Consider the first 
string equation 0 = V\L) n , n \ at order N~ 2 we have 


0 =|[iV 0 ] (s"±rd a d r + (s'f^rdldr + s'r'^d? 


r" ( Id 2 + i<9 r l + (r 1 ) 2 f 1 d 2 - 1 d r + ^d 2 d r 


6 


12 


+ [N 1 ]r / -d r + [N 2 ] f ■ [h°]V'{y/rh + s + y/rh x ) 


=U"^ 2 + K) 2 ^ 3 + U’z'^(j) 3 


+ z" ( ^0 2 + -=-Z“Vi ) + (-T ( ^ 


z'-z V02 + - 1 - Vi + U 2 (j)i + u\^ 


12 


J \2 




12 


1 _ 2 . 1 _! . 

—z Ei H-z U 3 

12 12 


(73) 


Proceeding similarly with the second string equation x = V'{L) n ^ n _ i we 
obtain 

0 =u" ( ~zcf) 2 + J + (u') 2 ( ^-zcj) 3 + -*-0 2 ) + u'z'\ip 3 


6 


12 


12 


12 


6 


6 


2 + (-') —03 - — Z *02 


12 


12 


+ V.\ ^“^01^ + U ' ^“«102^ + «201 + ^U?02 + -101 (74) 


We can solve for z\ in terms of lower genus quantities by taking the linear 
combination of equations zu'(73) + z'(74); this combination allows us to 
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isolate z\ using identity [65| . We obtain 


1 >< 2 , 1 _/ /n 3\ . , /. ty) 


=^3 ((*') + ^ («') j + 03 1 («') y + 12 


/\3 


+ 02 ( VX (V) 2 %' + ~Z ZU ' U " + o M l^ 


/„,// , 2 „/ 


12 

1 


6 


j\3 _/ _// 

+ 


12z 6 


+ 02 ( -zu"z' H—w?z?/ H— zu'z" -\ - u' (z') 2 

1 6 2 1 6 12 v ’ 

. u"z! u'z" u'(z') 2 1 . / 

1 12 12 12z 2 1 


(75) 


The appearantly complicated expression above simplifies in a fantastic way 


due to identities (65) and (66), which “unwind” the recurrence defining 0„ 


and 0 m . We would like to calculate which is related to z\ by F+ = 
d~ 2 {zi/z) — + log (z/x). It would therefore be nice if (75)/z could be grouped 


as a double exact derivative. Toward this purpose we recognize the following 


simplifications of blocks of terms that will appear in the equation (75 )/z\ 


/\3 


/ i , 1 , ^3 \ , A (y ,^2 /, yo 
*> 1 — + T2 (u) + *> i (M) z+ 1Tz 


,1 „ , 1 , „\ , (zf z’z" 1 , „ 

+ <h[- 6 uz + -uz ^ + i, 2 y-— + — + -uu 

=9. (yp'+y i 2 ' 


(76) 

(77) 


This simplification was achieved by applying identities (65) and (66) to the 


03 and 03 terms; collecting exact derivatives; then once again applying (65) 
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and ( 66 ). Similarly, 

02 


2 _/ 


U\Z 


+ 


(u') 2 z‘ 
2z ' 12 z 

+ (t)lUiu[ + -01 




u'z" 


1 2 , 0 ( 0) 2 
2^ v + ^r 

!( J\2 


+ 


u"z' u'(z') 


12z 12z 


12^ 2 


„ , 1 . 2 , U'z' 

= 4 (-*«, + *— 


Making these simplifications in (75); then using «i = ^u' and yet again 
identity ( 66 ), we obtain 

1 


Zl 


- =dx — 4>! u + —0 ! z + d x [ - fauf + 0 


12 


12 


■uV 


12 z 


- 01'Ul'Ul - 01 


z'u\ 


1 


=d x -—01 u" ~ ~(f>lZ + -01WJ + 01 


12 


12 


v! z' 


12 z 


2 z 
(78) 


=d x 0i 




A 2 


24 


u 


12 J ^ 12 

ft 


1 (u ') 2 0 — 2 z' z" + 2 zu'v! 

= 2 A dx (*') 2 - *(«') 2 

= - 4^ log [M 2 - Z («') 2 ] 

Finally we obtain the generating function for maps of genus 1: 

F w =ay^-0iogi 

z 12 x 

= 0 lo g [( Z 'f ~ Z{U') 2 } ~ 4 log*. 


(79) 


X 


(80) 
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3.0 Recursive calculation of z g and u g 

From lemma 3.2, it is clear that the highest genus terms in the string 
equations at order N~ 2g will be 

01 u 2g + z^ifiZg in equation 0 = [N~ 2a }V' (L) n ,„ 

0i« 2g + 0i z g in equation 0 = [N~ 29 ]V' {L)^^ 

Thus by the unwinding identity we can take combinations of equations to 
solve for the highest genus terms. 

zu'[N- 29 ]V\L) n)n + z'[N- 2a }V\L ) ntn _i isolates z g (81) 
z'[N- 2a }V’(L) n , n + 0[iV- 2 ^]l/'(L) nin _i isolates u 2g (82) 
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To calculate expressions for the functions u 2g+ \ in terms of lower genus quan¬ 
tities we do not need the string equations. It was shown in section 4.2.3 of 
[15] that 


2g+l R 

E -°m rvm 

7~T°® u 2g+l-m, 

m= 1 

and therefore the odd index functions u 2g+ \ are actually redundant in the 
sense that by applying the above identity several times u 2g+ \ can be ex¬ 
pressed purely in terms of lower genus u 2g with even index. (See [15]) for an 
explanation of the relation between u 2g and generating functions for maps 
with legs.) 

This proves that the functions z g and u 2g can be recursively calculated 
from the string equations. 

Definition 7. The polynomial degree of z or any x-derivative z^ is 1. The 
polynomial degree of u or any x-derivative u^ is The differential weight 
of any x-derivative u^ or z^ is m. 

Definition 8. A polynomial in u, z',u', z',u",... is homogeneous if all its 
terms have the same polynomial degree and the same differential weight. The 
polynomial degree of a product (or quotient) of homogeneous polynomials is 
given by adding (or subtracting) the polynomial degrees, and similarly for 
differential weights of products and quotients. 

Notice that the polynomial degree and differential weight are only de¬ 
fined for the subset of the rational functions in u, z' , u' , z', u",... satisfying a 
homogeneity property. Recall from unumiTH] the scaling relations 

z g (x, t ) =x 1 ~' 29 z g (1, (^ /2_ %)J =1 ) , (83) 

u g (x, t ) =x l/2 ~ 9 u g (l, (; x j/2 ~ 1 t j ) ( ] =1 ) . (84) 

A quantity with polynomial degree p and differential weight m has scaling 
structure x p ~ m G [fz^ 2 ~ l tj) d =1 ). Thus if one knows any two of the quantities 
polynomial degree, differential weight, and scaling structure, then the third 
is uniquely determined. 
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Theorem 2. From the string equations, we deduce that there exist formulas 
of the following forms for the functions z g and u g : 


u 2g 


U2g+1 


Pz g {u,z,u',z!,u",.. ■) 
((E) 2 - z{u') 2 f 9 ~* 
P U2g {u,z,u\z',u",...) 
{{z ') 2 - z{u'yf 9 ~ 3 

Pu 2g+1 (u,z,u\z',u",...) 
{{z'Y - z{u') 2 f 9 - 2 


(85) 

( 86 ) 

(87) 


where the functions P Zg , P U2g and P U2g+1 are polynomials in x-derivatives of u 
and z, homogeneous with respect to polynomial degree and differential weight. 
The differential weight of z g is 2 g, and the polynomial degree of z g is 1. The 
differential weight of u 2g is 2 g and the polynomial degree of u 2g is \. The 
differential weight of u 2g+ i is 2g + 1 and the polynomial degree of u 2g+ \ is f 


Let us comment that the above theorem statement is not sharp; consider 
for example the case of zp the formula we obtained has a denominator 
(( z') 2 — z(u') 2 ) 2 , but the theorem statement gives ((E ) 2 — z(u') 2 ) 5 . This does 
not contradict the theorem statement since the numerator and denominator 
are allowed to have common factors. This reason for the failure of sharpness 
is that each application of the unwinding identity cancels a common factor of 
( z') 2 — z(u') 2 from numerator and denominator. From the calculation of F b) 


in section 3.3| we see that the possibility of applying the unwinding identity 
is sensitive to the exact values of the string polynomials; since our formula 
for string polynomials is quite complicated, significant new ideas would be 
required to prove a sharp version of theorem [2] 

Proof of theorem [1| The case g — 1 is clear from the explicit calculation in 
section T3 Proceed by induction on g, assuming the theorem statement 
for all genus values less than g. First consider the right hand sides of the 
string equation 0 = [N~ 29 ]V'(L) n>n . Expressed in terms of the contributions 
of Motzkin paths, this is 


d 

{N- 2 ‘}V'(L) n:n = Y.(3 t i + ^) E < 88 ) 

j=1 
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where C(p) is defined in equation 32. From the scaling relations (83 84), it 
is clear that every term in this sum has scaling structure 


jV 2 - 2 »g ((j^-y)*). 


(89) 


We express equation (88) in the form guaranteed by theorem [lj then reduce 
all odd index functions U 2 h+i in terms of even index functions U 2 h /, and in 
each term apply one of the identities 


[N°h°]d™V'(h + s + rh- 1 ) =0 m+ 1 
[N^d^d^V’ih + s + rh- 1 ) =^Wi- 

The quantity [N-‘ 29 )V'(L) n:n is a sum of terms of the form 

(mf) (rn’,) im', ) (m 1 ’ ,) 

U2 9l ■ ■ ■ U 29k Z gk+1 . . . Z gK U 2g , . . . Z gWi . . . Z g , K/ 

X M (mi) . . . M (^) z (^+l) . . . z ^L) z r (0 m or ; 


(90) 

(91) 


where 


2 9 = ^2 2gi + 2g 'i + mi + m/ *' 


The highest genus terms are U 2 g (fti+z g z 1 ^i, but our inductive hypotheses 
apply to all other terms, which we will call forcing terms. Each forcing term 


has the scaling structure in equation (89), and by our inductive hypothesis, 
differential weight 2c/ — 1; also it is clear from the recurrence (64) that both 
(pm and have differential weight —1. Thus it follows that each term has 
polynomial degree One can see by the same argument that all forcing 
terms in the equation 0 = [N- 29 ]V'(L) n ^ n _i have polynomial degree zero and 
differential weight 2g — 1. By taking the linear combinations of equations 
given by (81), (82) we have that z g and u g are homogeneous with differential 


weight and polynomial degree as given in the theorem. 

It remains only to bound the power of (z') 2 — z{u') 2 appearing in denom¬ 
inators. From the equations at the end of section |2.4| we have 


P^{d s ,d r ) = Y^CtWt{d s ,d^ 


,)d l s W d l r {v) mod /, 


and the “mod /” cannot increase the total degree in d si d r . The sum is 
indexed by a tuple v, where the total |n| := Yhi v i i s l ess than or equal to 
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| A | + \rj\. From the proof of lemma |4j it is clear that the degree of the 
polynomial is at most |u| if this quantity is even, and |n| — 1 otherwise. 
Therefore the degree of P^(d s ,d r ) will be bounded by 


deg a.,e r P\°h( a *,dr) < |A| + |i)| + len(A) + len(jj). 


(92) 


Combining this estimate with lemma 3.2, we see that the index m of </> m or 
ijjm in (91) will be at most 


K + K' + L — 1 + ^ m, + m ■ 


Since </> m and V’m have denominator exponent 2m — 1, and are substituted 
into the string equations according to 


K =d™[h°]V\h + s + rh- 1 ), 

=d™~ 1 d^[h°}V\h + s + rh~ l ). 


and by our inductive hypothesis, the denominator exponent of forcing terms 


(91) will be at most 


2 + A + L — 1 + 'y ^ mi + m'^j — 1 + y ^ (&gi — 3 + 8 g[ — 3). 

This is clearly maximized by terms of either the form ( U 2 h or Zh)(u') p (z')^ p+ ^ m 
for which the bound is 8h — 4 + 4^mj; or terms of the form (u') p (z , )^ p+ ^ mi 
for which the bound is 8 h — 3. By the constraint (92), this is maximized 
in the 8h — 3 case. One can derive the same bound on the denominator 
exponents of forcing terms in the equation 0 = [N~ 29 ]V (L) nyn -i. Taking 
linear combinations of equations to isolate z g and U 2 g does not increase the 
denominator exponent. □ 


3.5. Formula for 


Following the procedure used in section 3.3 


one can compute F how¬ 
ever it is not feasible to perform this calculation by hand. Using a computer 
we have computed Z 2 by solving the string equations. We then computed 
F using the asymptotic Hirota-Szego relation 


a ( 2) =cu 2 




+ 



0 

x 


(93) 
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We find that the double integral d ~ 2 term can be computed in closed form. 
This is surprising, and we do not believe it to be a coincidence; however we 
are unable to prove that this will persist for higher values of g. We find that 


pi 2 ) 


+ 


240x2 ' 5760/ {z {u'f — (V) 2 ) 4 

- 24 (V) 10 + 96z {u'f {z'f + 24 zz" {z'f - 8// 3) {z'f - 144/ {u'f ( z'f 

- 84/ (-.u'f (z'f + 6z 2 {z'f {z'f - 96/ {u'f z" {z'f - 120 z 3 u'u™ {z'f 
+ 20// 4 ^ {z'f - 384/ {u'f u" {z'f + 384 z 3 u'u"z" {z'f 

- 84//V 3) {z'f + 172/ {u'f / 3) {z'f - 84/// 3) {z'f 

- 40/u'u (4) {z'f + 15/ {u'f {z'f + 64/ {z"f {z'f 

- 638/ {u'f {u"f {z'f - 340/ {u'f {z"f {z'f + 451z 3 {u'f z" {z'f 
+ 192/ {u"f z" {z'f + 48/ {ufu (3) {z'f + 252 z 4 u'z"u {3) {z'f 

+ 252///'/ 3) {z'f - 20/ {u'f / 4) {z'f - 256// (u") 3 {z'f 

- 768 z 4 u'u" {z"f {z'f + 1152/ {u'fu"z" {z'f - 168/ (/) W 3) {z'f 

- 152/ {u'f / 3) {z'f - 168/ {u'f z"z {3) {z'f + 80/ {u'f u (4) {z'f 

- 7/ (/) 8 {z'f + 384/ {u'f {z"f {z'f - 68/ {u'f {u"f {z'f 

- 430/ {u'f {z"f {z'f - 2/ {u'f z" {z'f + 1152/ {u'f {u"f z" {z'f 
+ 96/ (/)V 3) {z'f — 168/ {u'f z"u (3) {z'f — 168/ {u'fu"z^ {z'f 

- 20/ {u'f / 4) {z'f - 256/ {u'f {u"f z' - 768/ (/) 3 u" {z"f z' 

+ 252/ {u'f u"u {3) z' - 12/ {u'f z {3) z' + 252/ {u'f z" / 3) z! 

- 40/ {u'f u (A) z' + 64/ {u'f {z"f + 22/ {u'f {u"f - 4 / {u'f {z"f 
+ 7 / {u'f z" + 192/ {u'f {u"f z" - 24/ {u'f u {3) - 84/ {u'f z"u (3) 

- 84/ {u'f /'/ 3) + 20/ {u'f / 4) ^). (94) 


Appendix A. Identities for Motzkin paths 

In this section we prove some identities that were used above. The proofs 
treat r and s like formal parameters; that is, the lemma statements and 
proofs are equally valid with r and s replaced by u and z, respectively. 



Lemma 9 (Path reflecting, changing between monic and normalized OP’s). 
Let G be a smooth function. 


[h p ]G(h + s + rh x ) =r p ^ 2 [h p }G(y/rh + s + y/rh *) (A.l) 

=r p ^ 2 [h~ p ]G(y/rh + s + y/rh _1 ) (A.2) 

—r p [h~ p ]G(h + s + rh" 1 ) (A.3) 


Proof. An easy direct calculation. □ 

Lemma 10 (Path raising/lowering identity). Let G be a smooth function 
and p be an integer. Then we have the following identities, which correspond 
to lowering and raising the endpoint of a Motzkin path. 

d s [h p ]G(y/rh + s + y/rh" 1 ) —r p ^ 2 d r r 1 ^ 2 ~ p ^ 2 [h p ~ 1 ]G(y/rh + s + y/rh" 1 ) 

(A.4) 

d s [h p ]G(y/rh + s + y/rh" 1 ) —r~ p ^ 2 d r r p ^ 2+1 ^ 2 [h p+1 ]G(y/rh + s + y/rh" 1 ) 

(A.5) 


Iterating the above identities allows the right endpoint of a Motzkin path to 
be shifted to zero by 

d/ [h p ]G(y/rh + s + y/rh" 1 ) =/ p ^ 2 d/ [h°]G(y/rh + s + y/rh" 1 ) (A.6) 

Proof. First we prove ( |A.4 ) by direct calculation. 

[h p " 1 ]r p / 2 d r B 1 ~ p G(y/rh + s + y/rh" 1 ) =[h p ~ 1 ]r p ^ 2 d r G{h + s + rh" 1 ) 

= [h p " 1 }r p ^ 2 h" 1 G'{h + s + rh" 1 ) 

= [h p ]r p / 2 G' (h + s + rh" 1 ) 

= [h p ]G\y/rh + s + y/rh" 1 ) 

—d s [h p ]G(y/rh + s + y/rh" 1 ) 


The raising identity (A.5) now follows from the lowering identity (A.4) by 
applying [h r ]G(y/rh + s + y/rh" 1 ) = [h~ r ]G{y/rh + s + y/rh" 1 ) to both sides, 
then substituting p —» —p. □ 

Lemma 11 (Derivative swapping). Let G be a smooth function, then 

rd 2 [h°](y/rh + s + y/rh" 1 ) = ( d 2 - d r ) [h°](y/rh + s + y/rh" 1 ) (A.7) 
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Proof. 

dlfr[h°}G(y/rh + s + a / t / i -1 ) = [h°](h 2 + 2 + h~ 2 )G" (y/rh + s + a frh~ l ) 

=d 2 s ([ h ~ 2 ] + 2 [h°] + [h 2 ]) G{Vfh + s + yGh- 1 ) 
= ( rd 2 + 2d 2 + rd 2 ) [h°}G(^rh + s + a frhT 1 ) 


To get the last line above, we used the path raising/lowering identity ( A.6[ ). 
The lemma now follows by using d 2 ^ = 4 rd 2 + 2d r on the left hand. □ 


Lemma 12 (Integration by parts). 

a /rd s [h p ](h — h~ 1 )(yGh + s + a /rh~ l ) q =p[h p }(^rh + s + a /rh~ l ) q (A.8) 

Proof. Express the series coefficient as a contour integral about oo and inte¬ 
grate by parts. □ 
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Appendix B. Table of string polynomials 


A 

V 

p( a ) 


0 

0 

1 

0 

1 

0 

0 

-|rd r 

0 

1 

-d 

2 Ur 

0 

2 

0 

!rd d 

g ' L +s L/ r 

^rd 2 + T+r<9 r 

1 + 1 

0 

h rd l d r 

^rd 3 + ^rd s d r 

1 

1 

¥1 

|rd 2 d r 

0 

2 

+ d-Q 

6 u s ~ 12 r 

-rd d 

g' u s'-'r 

0 

1 + 1 

± r - 1 B 2 - +r _1 (9 + —d 2 d 

12 U s 12 Ur ' 12 u s u r 

d-d 3 — +d d 

3 

0 

0 

-+rd 2 

12 u s 

2 + 1 

0 

0 

1 rd 3 

6 s 

1 + 1 + 1 

0 

0 

— —rd 4 

24 1 s 

2 

1 

—(9 3 

12 u s 

-d- r d 2 d 

12 u s u r 

1 + 1 

1 

d-d 4 

24 s 

—A rd 3 d 

12 U s u r 

1 

2 

—(9 3 

12 

-d- r d 2 d 

12 ' u s u r 

1 

1 + 1 

—d 3 d 

12 u s u r 

-d-d 4 + +d 2 d 

24 ~ 24 u s u r 

0 

3 

d-d 2 

12 w s 

0 

0 

2 + 1 


0 

0 

1 + 1 + 1 

~ 14'' 'I's'l 

0 
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